SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

FERROELECTRICS AND 
MULTIFERROICS 

MAGNETIC PROPERTIES AND 
MATERIALS 



Received 
3 October 2013 

Accepted 
20 January 2014 

Published 
5 February 2014 



Correspondence and 
requests for materials 
should be addressed to 
J.-M.L. (Iiujm@nju.edu. 

cn) 



Experimental observation of 
ferrielectricity in multiferroic DyMn205 

Z. Y. Zhao 1 , M. F. Liu 1 , X. Li 1 , L. Lin 1 , Z. B. Yan\ S. Dong 2 & J. -M. Liu 1 

1 Laboratory of Solid State Microstructures, Nanjing University, Nanjing 2 1 0093, China, 2 Department of Physics, Southeast 
University, Nanjing 2101 89, China. 

One of the major breakthroughs associated with multiferroicity in recent years is the discovery of 
ferroelectricity generated by specific magnetic structures in some magnetic insulating oxides such as 
rare-earth manganites RMn0 3 and RMn 2 0 5 . An unresolved issue is the small electric polarization. 
Relatively large electric polarization and strong magnetoelectric coupling have been found in those 
manganites of double magnetic ions: magnetic rare-earth R ion and Mn ion, due to the strong R-Mn (4f-3d) 
interactions. DyMn 2 0 5 is a representative example. We unveil in this work the ferrielectric nature of 
DyMn 2 0 5 , in which the two ferroelectric sublattices with opposite electric polarizations constitute the 
ferrielectric state. One sublattice has its polarization generated by the symmetric exchange striction from the 
Mn-Mn interactions, while the polarization of the other sublattice is attributed to the symmetric exchange 
striction from the Dy-Mn interactions. We present detailed measurements on the electric polarization as a 
function of temperature, magnetic field, and measuring paths. The present experiments may be helpful for 
clarifying the puzzling issues on the multiferroicity in DyMn 2 0 5 and other RMn 2 0 5 multiferroics. 

Multiferroics have been intensively investigated for ten years since the pioneer works on BiFe0 3 1 and 
TbMn0 3 in 2003 2 . In particular, the discovery of magnetically induced ferroelectrics (the so-called type- 
II multiferroics) has comprehended our understanding of multiferroicity 2 4 . In these materials, electric 
polarization P is believed to be correlated with particular magnetic orderings below certain temperatures and thus 
the cross-coupling between ferroelectricity and magnetism is significant, allowing possible magnetic control of 
ferroelectricity or/ and electric control of magnetism 511 . To dates, what keeps the research interest alive is the 
possibility of unveiling microscopic physics which is substantially different from our earlier knowledge and even 
general principles for guiding the design and synthesis of multiferroics of promising practical applications 51215 . It 
is noted that most discovered type-II multiferroics so far either have low ferroelectric and magnetic transition 
temperatures or exhibit small electric polarization/weak magnetization. 

While conventional ferroelectrics exhibit the electric polarization via the structural symmetry-breaking transi- 
tions from high symmetric paraelectric (PE) phase 16 , for those type-II multiferroics the primary order parameter 
is magnetic rather than structural. Two major magnetic mechanisms for the ferroelectricity generation have been 
proposed. One is the asymmetric exchange striction, in which the inverse Dzyaloshinskii-Moriya (DM) inter- 
action associated with the non-collinear spin ordering drives the structural symmetry-breaking 617,18 . The other is 
the symmetric exchange striction, in which specific collinear spin ordering, such as the E-type antiferromagnetic 
(AFM) ordering 19,20 and ttl-l ordering 21 , drives the structural symmetry-breaking. 

Interestingly, another class of multiferroics, in which the symmetric exchange striction is believed to play major 
roles, is rare-earth manganites RMn 2 0 5 family 22 . The RMn 2 0 5 family shows complicated lattice distortions and 
spin structures, and exhibits multifold competing interactions, large electric polarization, and remarkable mag- 
netoelectric (ME) responses 5-16,23 . Nevertheless, partially due to the multifold competing interactions, the multi- 
ferroic transitions and underlying mechanisms in RMn 2 0 5 are not yet well understood 8,9 . All members of this 
RMn 2 0 5 family have similar structural ingredients 16 . The lattice structure projected on the afo-plane is shown in 
Fig. 1. The Mn ions are partitioned into Mn 3+ and Mn 4+ , which are coordinated respectively in square pyramid 
Mn-0 units and octahedral Mn-O units. On the afo-plane, the octahedra and pyramids are corner-sharing by 
either the pyramid bases or pyramid apex, and the adjacent pyramids are connected with their bases. Along the 
c-axis, the octahedral sharing edges constitute linear chains. Each Mn 3+ ion is located in between two Mn 4+ ions, 
and the R 3+ ions are located on the alternative layers between two Mn 4+ ions. 

In RMn 2 0 5 , the Mn spin interactions are characterized by the three dominant components J 3 , J 4 , and J 5 , plus 
additional long-range components 16,24 . Their competitions lead to consecutive commensurate antiferromagnetic 
(C-AFM) and incommensurate AFM (IC-AFM) ordering sequence 25,26 . Also, the 3d-4/ interactions can't be 
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Figure 1 | A schematic drawing of the lattice structure of DyMn 2 0 5 with 
the three major Mn-Mn spin interactions J 3 , ] 4 , and J 5 . The ions and 
coordinates are drawn for guide of eyes. 

neglected if the R ion has big moment, and the strong R-Mn coupling 
allows even more fascinating spin structure evolution. The ferroelec- 
tricity and its magnetic origins are thus far from understanding. 
Here, we choose DyMn 2 0 5 as a representative example for illustra- 
tion 8,9,24 . The paramagnetic phase above temperature T ~ 43 K tran- 
sits into an IC- AFM phase, followed by a C- AFM phase below T Nl = 
—40 K, and then by the coexistence of an IC-AFM phase and a C- 
AFM phase below T N2 ~ 28 K. This coexistence is again replaced by 
two coexisting IC-AFM phases below T N3 — 20 K. At T < T Dy — 8 K, 
the Dy 3 * spins order independently. The structural and interaction 
origins for these magnetic transitions were discussed extensively, 
while no full consistency has been reached 8,9,25 . 

With respect to the magnetic structures, our understanding of the 
ferroelectricity is even in the earlier stage. While it is believed that 
the C-AFM phase is ferroelectric and the IC-AFM phase is not, the 
measured results are not always consistent with this prediction 8,9 . 
Basically, the measured P most likely aligns along the fo-axis. 
However, the measured data in experiment by Hur et al show that 
the P appears below T N1 (T FE1 ) and changes its sign from negative 
value to positive one at a certain T lower than T N2 ~ 27 K, giving a 
feature of ferrielectric (FI) state 8 . However, in experiment by 
Higashiyama et aP, the measured P(T) experiences several transi- 
tions which correspond one-to-one to the magnetic transitions, and 
the system becomes non-ferroelectric below T Dy ~ 8 K (the so-called 
X-phase) while the ferroelectric nature of this X-phase remains 
unclear. Besides, the magnetic origins for these transitions were dis- 
cussed in details 16 . So far available data on the ME effect of RMn 2 0 5 
family are also materials-dependent. Remarkable ME response was 
observed in the low- T range for those materials with big 4/ magnetic 
moments, where the 3d-4f (R-Mn) coupling is strong enough in 
determining the magnetic structure 8,9,12,26,27 . 

The inconsistencies and insufficient data, as partially highlighted 
above, suggest a critical appealing for revisiting the electric polariza- 
tion and its response to magnetic field in RMn 2 0 5 (here DyMn 2 0 5 ). 
On the basis that DyMn 2 0 5 has strong Dy-Mn interactions in addi- 
tion to the dominant Mn-Mn interactions 28 , one has reasons to 
expect a ferrielectric state with more than one polarization compon- 
ent. With no doubt, convincing evidences with this ferrielectric state 
become primarily critical for understanding the multiferroicity of 
DyMn 2 0 5 and more generally the RMn 2 0 5 family. In this work, it 
will be suggested that DyMn 2 0 5 is a ferrielectric composed of two 
anti-parallel ferroelectric sublattices. The electric polarizations of the 
two sublattices have different microscopic origins, with one arising 



from the Mn-Mn symmetric exchange striction and the other from 
the Dy-Mn symmetric exchange striction. It is noted that the ME 
effect can be reasonably explained by this ferrielectric model. We 
employ a modified pyroelectric current (mPyro) method to track 
the evolution of the electric polarization upon various paths, while 
a detailed discussion on the methodology for measuring the polar- 
ization and a description of this mPyro method can be found in the 
Supplementary document. 

Results 

Multiferroic phase transitions. We first look at the phase transition 
sequence in terms of specific heat Cp, magnetization M, and dielectric 
constant 6 as a function of T, as shown in Fig. 2. For reference, the 
released current I tot {T) (i.e. Ip yr0 (T)) using the mPyro method at a 
warming rate of 2 K/min is shown in Fig. 2(d), while a demonstra- 
tion of the mPyro method in precisely measuring the pyroelectric 
current released from polarized charges is given in the Supplemen- 
tary document. The Cp(T) curve shows clear anomalies roughly at 
T N1 ~ 40 K, T N2 ~ 27 K, and T Dy ~ 8 K, while the peak at T N3 ~ 
20 K, if any, is weak. These anomalies reflect the sequent magnetic 
transitions from the paramagnetic phase to the C-AFM phase, to the 
coexisting IC-AFM phase plus C-AFM phase, then to the two IC- 
AFM coexisting phases, and eventually to the independent Dy 3 * spin 
order plus IC-AFM phase, consistent with earlier reports 26 . However, 
no features corresponding to these transitions, except the indepen- 
dent Dy 3 * spin ordering at T Dy , were observed in the measured M-T 
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Figure 2 | (a) Specific heat normalized by temperature (C P IT), (b) 
magnetizations (M) under the ZFC and FC conditions, (c) dielectric 
constant (s), and (d) released current (/,„, = I pyTO ) by the mPyro method 
at a warming rate of 2 K/min with a poling electric field E poIe =10 kV/cm, 
as a function of T, respectively. The dielectric constant was measured at 
frequency of 100 kHz with a bias of 50 mV, and no remarkable frequency 
dispersion was observed. The phase regions proposed in literature are 
labeled on the top. 
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data, mainly due to the fact that the Dy 3 " 1 " moment is much bigger 
than the Mn 3+ /Mn 4+ moments. The anomalies in the s-T curve at 
these phase transitions reflect the magneto-dielectric response, as 
revealed earlier 8 . Interestingly, a series of anomalies in the I pyro -T 
curve at these transition points are available, as shown in Fig. 2(d), 
consistent with earlier report 8 too, evidencing the strong ME effect. 
In addition, these magnetic transitions may be path-dependent, and 
the features in the cooling run are different from those in the 
warming run. Our data on the electric polarization below also 
illustrate this dependence. 

Nonzero polarization of the X-phase. The measured I pyr0 -T curves 
and as-evaluated P-T curves under E po i e = 10 kV/cm are plotted in 
Fig. 3, given different starting temperatures (T e „ d ), where the T end is 
the temperature to which the sample is cooled down from high-T 
paramagnetic state under E po \ e =10 kV/cm. It is seen that the Ip yro - T 
curves exhibit clear anomalies at the magnetic transition points (T N1 , 
T N2 , T N3 , and T Dy ), and the P- T curve at T end = 2 K is similar in shape 
to that reported in Ref. 8. 

We suggest that the X-phase is ferrielectric (i.e. ferroelectric in 
general sense). First, the measured data in Fig. 3(a) and (b) indicate 
that the X-phase has nonzero electric polarization. The pyroelectric 
current in both the X-phase region and the other three ferroelectric 
regions is much bigger than 0.3 pA, the background level in the 
present experiment. It is also observed that the pyroelectric current 
depends on the poling field, the bigger the current the higher the 
field. The mPyro method used in the present experiment is different 
from the Pole method used in Ref. 9 for the electric polarization 
measurement (see the Supplementary document). In the Pole 
method, the current flowing across the sample under a relatively 
low electric field is measured during the sample cooling. In this case, 
by assuming that the leakage current at low Tis much lower than the 
polarization current, one may evaluate the electric polarization from 



the polarization current data directly. As seen in Ref. 9, the measured 
polarization in the X-phase region was indeed negligible, suggesting 
that the X-phase is non-ferroelectric. Clearly, if the X-phase is anti- 
ferroelectric, or ferrielectric with two comparable antiparallel polar- 
ization components, the polarization current can be small and even 
comparable with the leakage current. In this sense, the Pole method 
may not be applicable for identifying the ferrielectric or antiferro- 
electric state. 

Second, each measured P(T) curve shown in Fig. 3(a) ~ (d) indi- 
cates a negative-positive sign change with decreasing T, suggesting 
immediately that DyMn 2 0 5 is a ferrielectric (FI) rather than a nor- 
mal ferroelectric 8 . The sign change would be the consequence of 
competition between the two ferroelectric sublattices whose polar- 
ization components should exhibit different T-dependences. In this 
case, careful measurement using the mPyro method can provide 
critical data on details of the ferrielectric state and the different In- 
dependences. 

Path-dependent polarization. Given a fixed E po i e , the I pyr0 (T) and 
thus the P(T) show remarkable T end - dependent behaviors, i.e. path- 
dependent. To see clearly this path-dependence, the I pyro -T curve 
with T ent i = 2 K is shown by a thin dashed line in each plot of 
Fig. 3. The temperature T P = 0 , at which the P(T) changes its sign, 
is plotted as a function of T end in Fig. 4(a). The T P = 0 shifts ~9 K 
when T end increases for ~ 10 K, implying that the ferrielectric state is 
not robust against thermal fluctuations (T) or external field (£ po / e ). 
Supposing that the ferrielectric state is composed of two sublattices, 
one expects that at least one of them is strongly T-dependent or £ po / e - 
dependent. What surprising us is the Ip yr0 -T curves as T end > 12 K, 
some of which are plotted in the right column of Fig. 3. In spite of 
positive Ep 0 i e , the measured P data below T N1 are negative. The 
negative I pyro peak at T N1 remains nearly unchanged even with 
Tend = 38 K, very close to T N1 = 40 K. Such a negative P can't be 
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Figure 3 | Measured pyroelectric current I pyro = I P and evaluated electric polarization Pas a function of Tat T end = 2 K (a), 6 K (b), 8 K (c), 10 K (d), 
15 K (e), 25 K (f), 33 K (g), and 38 K (h), respectively. The warming rate is 2 K/min. For reference, the ip-T data at T end = 2 K are inserted. 
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Figure 4 | (a) Evaluated crossing temperature T P = 0 at which the measured 
P{ T) changes its sign, as a function of T end . (b) Evaluated P{ T) curve and 
P( T e „ d ) curve. The warming rate for the pyroelectric current probing is 
2 K/min. 

possible in a normal ferroelectric, unless the P has two components 
which are anti-parallel to each other. 

The T e „ d -dependences of the Ip yr o(T) and P(T) curves are strange 
at the first glance. It is reflected that the electric polarization has the 
magnetic origin, since no structural phase transitions occur below 
T w .We present in Fig. 4(b) the measured P value at T end , i.e. P(T enc i), 
where the P(T) curve with T md = 2 K is inserted for comparison. For 
a normal ferroelectric, the P{T end ) should overlap with the P(T) by 
setting T end = T. Here, the overlapping only occurs below T Dy , noting 
that the P(T end ) is always larger than the P(T). The difference 
between them maximizes at T end = T ~ 10 K and —24 K and 
becomes negligible as T end —> T N1 , suggesting that the magnetic 



transitions below T N1 have the path-dependent characteristic, while 
the first-order or second-order nature of these transitions deserves 
for additional clarification. In fact, combining the P{T end ) and P(T) 
data generates a double-loop like hysteresis, as shown in Fig. 4(b). 

Discussion 

For DyMn 2 0 5 , the symmetric exchange striction effect arising from 
the specific Mn/Dy spin alignment is the main mechanism for the 
electric polarization 16 . For a simplification consideration, we don't 
take into account the contributions from the noncollinear spin 
orders. However, the effect of the independent Dy spin ordering at 
T Dy imposes significant effect on the electric polarization, due to the 
strong Dy-Mn interactions, and thus will be considered. 

Referring to relevant literature on DyMn 2 0 5 24 , we present in 
Fig. 5(a) the spin structure projected on the afo-plane over the T- 
range between T N1 and T Dy . The square pyramidal and octahedral 
structural units surrounding the Dy 3+ , Mn 3+ , and Mn 4+ spins are 
drawn for a better view in Fig. 5(b). The light gray and gray structural 
units shift 1/4 lattice unit from each other along the c-axis 16 . Along 
the fo-axis, one finds two types of three-spin blocks each centered on a 
Mn 4+ spin, as shown in Fig. 5(c) and (d), respectively. One is block A, 
consisting of one Mn 4+ -0 octahedron connected with two pyramid 
units each with one Mn 3+ spin inside (Fig. 5(c)). The other is block B, 
consisting of one Mn 4+ -0 octahedron connected with two Dy 3 * 
spins located in the space surrounded by the Mn0 6 and Mn0 5 units 
(Fig. 5(d)). Because of the symmetric exchange striction, the two 
Mn 3+ ions in the block A shift roughly up and the two Dy 3 * ions 
in the block B shift down with respect to the Mn 4+ ions. Therefore, 
one electric polarization component (Pmm) m the block A and one 
polarization component (Pdm) in the block B are generated. They are 
roughly anti-parallel to each other but align along the fo-axis. The 
whole lattice as the consequence of the alternating stacking of the two 
types of blocks is therefore a ferrielectric lattice composed of two FE 
sublattices. 

Different from the P MMl the P DM originates from the Dy-Mn 
interactions and thus depends on the Dy 3 * spin order. Above T Dy , 
the Dy 3 * spins may order in coherence with the Mn spin ordering 
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Figure 5 | Proposed spin structure at a temperature lower than T N2 and higher than T Dy , referring to neutron scattering data available in literature, (a) 

The spin structure projected on the ub-plane with the square pyramidal Mn 3+ -0 2 ~ unit and octahedral Mn 4+ -0 2 ~ unit shown in (b). The structural block 
A, composed of one Mn 4+ -0 2 ~ octahedra connected by two Mn 3+ -0 2 ~ pyramids roughly along the fo-axis, is shown in (c). The structural block B, 
composed of one Mn 4+ -0 2 ~ octahedra connected by two Dy 1+ roughly along the fo-axis, is shown in (d). The proposed polarizations P MM and P DM 
generated by the two types of blocks due to the symmetric exchange strictions, are labeled in (c) and (d), respectively. 
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around T N1 or T N2 , due to the strong Dy-Mn interactions. At T < 
T Dy , this induced Dy 3+ spin ordering is partially and gradually 
replaced by the independent Dy 3 * spin ordering, although the details 
of the Dy 3+ spin ordering sequence has not been well understood. To 
this stage, one has reason to argue that the P DM will show much more 
significant T-dependence than the Pmm- The reason can be discussed 
considering the weak ordering of the Dy 3 * spins themselves. This 
ordering is sensitive to the 3d-4f interactions and external field. We 
discuss this issue below from various aspects. 

First, the 4/ interactions in some transition metal oxides are quite 
localized and the R 3+ spins alone can't order unless the temperature 
is very low (the ordering point T R is less than ~2 K), as identified in 
oxides without 3d moments but only the 4/ moments, as seen in 
R 2 Ti 2 0 7 etc with R = Gd, Tb, Ho, and Er etc 29 , where the Ti 4+ has 
no magnetic moment and thus no 4f-3d interaction is available. In 
some other oxides, the R 3+ spins can't order even at extremely low 
temperature, leading to spin liquid or spin ice states due to the crystal 
fields and quantum fluctuations 30,31 . In these oxides, a magnetic field 
of —l.OTesla is sufficient to break the original spin orders and 
enforce the parallel spin alignment. 

Second, for DyMn 2 0 5 and other RMn 2 0 5 /RMn03 with R = Gd, 
Ho, Er etc, the situation can be different since the R 3+ spins coexist 
with the Mn spins. The independent R 3+ spin ordering can occur at a 
T R as high as 6 ~ 10 K (here T Dy ~ 8 K). Furthermore, if the Mn spin 
ordering occurs well above T R , an additional R 3+ spin ordering at a 
temperature higher than T R will be induced by the Mn spin orders 
due to the 4f-3d interactions. These observations suggest that the 
4f-3d interactions can enhance the T R value if the sign of interaction 
is consistent with that of the 4/-4/ exchange interaction 16,32,33 . A typ- 
ical case is seen in DyMn0 3 , where the Mn 3+ spins order antiferro- 
magnetically at —38 K and are locked in the noncollinear spiral order 
at ~20 K. Slightly below this locking point, the Dy 3 * spins order in a 
coherent manner with the Mn 3+ spin order. This induced Dy 3+ spins 
order sustains until T Dy ~ 7 K at which the independent Dy 3 * spin 
ordering enters 34 . In addition, a magnetic field of 1.0 ~ 2.0Tesla is 
sufficient to break the coherent and independent Dy 3 * spin orders, 
while much higher field is needed to melt the Mn spin orders. 

The above discussion suggests that the Dy 3+ and other rare-earth 
moments have relatively weak 4/-4/ exchange-coupling with respect 
to the 3d moments such as Mn spins here. This discussion thus serves 
as the model basis on which the Pom an d Pmm as a function of T are 
evaluated, respectively. For simplification, the effect of independent 
Dy 3+ spin ordering below T Dy on the Mn spin order is assumed to be 
weak if any. It can be reasonably assumed that the Pmm initiating at 
T N1 increases rapidly in magnitude with decreasing T and becomes 
saturated in the low T range, because the Mn spin order is already 
well developed below T N1 . Consequently, the P DM as a function of T 
can be extracted. Take the data with T end = 8 K ~ T Dy as an example. 
The measured P(T) data are plotted in Fig. 6(a). The Pmm(T) curve is 
extracted based on the above assumption, and then Pom(T) = P(T) 
— Pmm(T) is evaluated. For a clear illustration, the two ferroelectric 
sublattices on the afo-plane are schematically drawn in Fig. 6(c) and 
(d), and a combination of them constitutes the ferrielectric lattice in 
Fig. 6(b). As expected, the P DM increases gradually with decreasing T 
until T ~ 20 K, below which a much more significant T-dependence 
than that for the Pmm is then exhibited. 

As Tend = 2 K < T Dy , the effect of the independent Dy 3 * spin 
ordering on the P DM takes effect. At T = T en ^ some Dy 3+ spins are on 
the track of the independent ordering, leading to disappearance of 
Pom at some lattice sites. The Pom sublattice is thus partially melted 
away, giving rise to a smaller Pom- This is the reason for the low T P = 0 
and small \P\ below T P = 0 , with respect to the case of T enc i = 8 K. 
Here it should be mentioned that the difference in the Pom(T) curve 
between the case of T eni j = 8 K and that of T ent i = 2 K reflects the 
difference in the magnetic structures between the two cases. The origin 
lies in the fact that the magnetic transition at Toy is path-dependent. 



0.6 
0.4 
0.2 

0.0 



(a) 

P 


v P 

\ DM 


■ / P \^\_ 

/ MM ^T^V 







10 



a 



20 30 
T(K) 



40 



50 



'l St 





;• 




•t— «t- 


'i 

;« 




•t- -K 


if 





p IP tp 

MM DM 

(b) (c) (d) 

Figure 6 | (a) Evaluated electric polarizations P MM and P DM from the two 
ferroelectric sublattices of the proposed ferrielectric model, as a function of 
T, where P = P DM + P MM and T enl j = 8 K. The proposed ferrielectric lattice 
and the associated two sublattices, all projected on the afo-plane, are 
schematically drawn in (b), (c), and (d), respectively. 

Given the ferrielectric model shown in Fig. 5 and the different 
Pdm(T) and Pmm(T) behaviors, a puzzling issue appears: why does 
the measured Pmm (° r P) remain negative even though T end is higher 
than T N2 ? For T ent i > T N2 , the Pom should be much smaller than the 
Pmm an d thus a poling by a positive E po i e would generate a positive 
Pmm- In this case, the measured I pyw and P should be positive, contra- 
dicting with the measured data. At this stage, we have no convincing 
explanation of this anomalous phenomenon. One possible reason is 
that the Pom(T) is sensitive to the E po \ e . Considering the fact that the 
Dy 3 * spins have weak exchange coupling, as addressed above, one 
expects that the electric field driven alignment of the Dy 3+ spins 
coherently with the Mn spins would be energetically easy. 
Therefore, the Pom can De remarkably enhanced by the £„ 0 ; e . If it 
is the case, the electric poling during the cooling sequence can 
enhance the Pom remarkably while the Pmm is roughly unchanged, 
so that the Pom around T end is larger than the Pmm in magnitude. 
This results in the alignment of the Pmm opposite to both the fL ofe 
and Pom- After the removal of the E po i e at T en ^ the Pom shrinks back 
to a value smaller than the Pmm- Consequently, the pyroelectric 
current remains negative. Another possible explanation for this 
strange phenomenon is the ferroelastic effect in DyMn 2 0 5 35 , which 
makes the Pmm domains be clamped along a direction opposite to the 
E po i e during the poling process. However, this assumption remains to 
be confirmed. 

Obviously, referring to the ferrielectric model, one immediately 
predicts that the ME parameter AP(fl) = [P(H) - P(H = 0)] is 
negative and also remarkably T-dependent, since a magnetic field 
as big as 1.0-2.0Tesla is sufficient to align the R 3+ spins even at an 
extremely low T 9 . For DyMn 2 0 5 here, this effect suggests that a 
magnetic field of 1.0-2.0 Tesla re-aligns the Dy 3 * spins along the 
field direction, while the Mn spins remain robust. Therefore, the f fj 
or HI pattern in the block B is broken, as shown in Fig. 7, leading to 
the P D m ~ 0- The spin structure in the block A remains roughly 
unchanged and thus does the Pmm- The experimental data conform 
this prediction, as presented in Fig. 8. In Fig. 8(a) are plotted the P-T 
data at P end = 2 K, where the Pom(T) and Pmm(T) under H = 0 are 
presented too. The AP{H > 2 T) should not be much less than the 
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Figure 7 | Schematic drawing of the spin alignments in the block A and 
blockB, respectively, under a downward magnetic field Hat T< T Dr The 
Dy" + spins can be easily re-aligned by H while the Mn spins can't, implying 
that the P DM = 0 at T < T Dy . 

Pdm in magnitude although their signs are opposite. Our data also 
support this prediction. The ferrielectric state as a basis for this ME 
effect is then confirmed. Here, it should be mentioned that the model 
shown in Fig. 7 assumes that magnetic field H is parallel to the fc-axis. 
Nevertheless, our samples are polycrystalline and the grains/mag- 
netic domains are randomly oriented. Taking into account of the 
polycrystalline nature, one still can expect that the |f J, and HI 
patterns of most block B units in the sample will be broken. This 
model explanation is thus qualitatively reasonable. 

To this stage, we have presented a qualitative explanation of the 
major features associated with the electric polarization and ME 
effect in DyMn 2 0 5 , based on the proposed ferrielectric model. 
Nevertheless, several issues remain yet unclear or unsolved: (1) No 



detailed discussion on the possible ferroelectric phase transitions at 
the magnetic transition points T N2 , T N3 , and even T Dy , respectively, 
has been given. (2) The path-dependence of the electric polarization 
is attributed to the magnetic transitions which are path-dependent. 
The first-order ot second-order nature of these magnetic transitions 
remains to be clarified. (3) An uncertain point regarding the present 
ferrielectric model is the response of the Pdm to electric field which is 
assumed to be remarkable in order to account for the experimental 
observations. Searching for convincing evidence on this assumption 
is challenging although the assumption itself is physically reasonable. 
A careful characterization of the Dy 3 * spin structures at various T is 
critical for dealing with these issues. It was reported that the element 
specific X-ray resonant magnetic scattering (XRMS) is a powerful 
tool although neutron scattering may face the problem of large 
absorption by Dy nuclei 36 . 

In summary, extensive multiferroic measurements on DyMn 2 0 5 
have been carried out, and the complicated electric polarization 
behaviors have been characterized. It is revealed that the electric 
polarization in DyMn 2 0 5 does consist of two antiparallel compo- 
nents, demonstrating the ferrielectric state at low temperature. The 
two electric polarization components are believed to originate from 
the symmetric exchange striction. One is generated from the 
Mn 3+ -Mn 4+ -Mn 3+ blocks with the J. T T and Til spin alignments, 
which is robust against temperature and magnetic field. The other 
is generated from the Dy 3+ -Mn 3+ -Dy 3+ blocks with the J. j.f and ffj. 
spin alignments, which is sensitive to temperature and magnetic 
field. The present work represents a substantial step towards a full- 
scale understanding of the electric polarization in DyMn 2 0 5 and 
probably other RMn 2 0 5 family members. 

Methods 

Polycrystalline DyMn 2 0 5 samples were used for the present experiments. The sam- 
ples were prepared by standard solid state sintering. Stoichiometric amount of 
Dy 2 0 3 (99.99%) and Mn 2 0 3 (99%) was thoroughly mixed, compressed into pellets, 
and sintered at 1200 L, C for 24 h in an oxygen atmosphere with several cycles of 
intermediate grindings. For every sintering cycle, the samples were cooled down to 
room temperature at 100 'C per hour. The as-prepared samples were cut into various 
shapes for subsequent microstructural and property characterizations. The sample 
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Figure 8 | Measured ME responses and proposed model, (a) The measured P(T) curves and proposed Pmm(T) and Pdm(T) curves under H = 0 and 
ff 3> 0 (e.g. ~2 T). It is suggested that the P MM is robust against H while the P DM can be seriously suppressed by H, due to the field induced Dy 3+ spin 
realignment as proposed in Fig. 7. The ferrielectric lattice at H = 0 is shown in (b), which transfers into the lattice in (c) at H 2> 0. This lattice in (c) is 
composed of the P MM sublattice shown in (d) plus the P DM sublattice shown in (e). P = P DM + P MM . 
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crystallinity was checked using X-ray diffraction (XRD) with Cu Ka radiation at room 
temperature. 

Measurements on the specific heat (Cp), magnetization (M) and dc magnetic 
susceptibility (x), dielectric susceptibility (e) and electric polarization (P) of the 
samples were carried out. All the data presented in this work were obtained from the 
polycrystalline samples. The M and / were measured using the Quantum Design 
Superconducting Quantum Interference Device (SQUID) in the zero-field cooled 
(ZFC) mode and field-cooling (FC) mode, respectively. The cooling field and mea- 
suring field are both 1000 Oe. The C p was measured using the Quantum Design 
Physical Properties Measurement System (PPMS) in the standard procedure. 

The polarization P was measured using the modified pyroelectric current (mPyro) 
method with different starting temperature T enc t — 2 K — 38 K, respectively. Each 
sample was polished into a thin disk of 0.2 mm in thickness and 10 mm in in-plane 
dimension, and then sandwich -coated with Au layers as top and bottom electrodes. 
The measurement was performed using the Keithley 6514 A and 6517 electrometers 
connected to the PPMS. In details, each sample was submitted to the PPMS and 
cooled down to —100 K. Then a poling field E po \ e —10 kV/cm was applied to the 
sample until the sample was further cooled down to T enc j, at which the sample was 
then short-circuited for sufficient time (>30 min) in order to release any charges 
accumulated on the sample surfaces or inside the sample. The recorded background 
current noise amplitude was —0.3 pA. Then the sample was heated slowly at a 
warming rate up to a given temperature T 0 — 60 K > T N , during which the released 
current I tot was collected. Similar measurements were performed with different 
warming rates from 1 K/min to 6 K/min and the collected I tot data are compared to 
insure no contribution other than pyroelectric current I pyro . Finally, polarization P( T) 
was obtained by integrating the collected I pyro {T) data from T 0 down to T m ^, The 
validity of this procedure was confirmed repeatedly in earlier works' 5 and the data 
presented in the Supplementary document. 

In addition, the e data at various frequencies as a function of T were collected using 
the HP4294A impedance analyzer with an ac-bias field of —50 mV. Besides the e-T 
data and P- T and data, we also measured the response of P to magnetic field H in two 
modes. One is the isothermal mode with which the variation in P in response to the 
scanning of H was detected and the other is the iso-field mode with which the P-Tdata 
under a fixed H were collected. By such measurements, one can evaluate the ME 
coupling of the samples. We define AP{H) — P{H) — P{H — 0) as the ME parameter. 

We also employed the PUND method to obtain the P-E loops at various tem- 
peratures, using the identical procedure as reported in literature e.g. Ref. 37. Our data 
are quite similar to reported ones from other groups (see the Supplementary docu- 
ment). 
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